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1. INTRODUCTION 

In the field of electrical energy, there are many different and varied sources, where we find 
renewable and non-renewable sources. Wind energy is a type of renewable energy. The latter has become the 
main pillar of the global economy and the production of electrical energy from the wind energy we use wind 
turbines. Currently, there are two types of turbines to generate electric power, which are single-rotor wind 
turbine (SRWT) and dual-rotor wind power (DRWP). But, the SRWT is widely used in electric power 
generation compared to DRWT and this is due to that DRWT is a newly discovered compared to the SRWT 
system. Also, the generation of electrical energy depends on converting mechanical energy into electrical 
energy, using generators, as there are several generators, for exemple: synchronous generator (SG), 
permanent magnet synchronous generator (PMSG), and doubly-fed induction generator (DFIG). In our 
article, we will rely on generating electric energy using a DFIG-based DRWT system. In the field of 
scientific resarche, there are several methods for controlling generators, for example: indirect vector control 
(IVC) [1]-[3], backstepping command [4]-[6], direct torque command (DTC) [7]-[10], direct vector 
command [11]-[13], direct power command (DPC) [14]-[18], intelligent command [19]-[22], non-linear 
command [23]-[25] and hybrid command [26]-[30]. 

In this work, the DPC technique with the application of the genetic algorithm (GA) and traditional 
PWM technique has been considered. The original contribution of this work is the application of the GA 
method in the DPC system with three-phase DFIG-based DRWT and simulation investigation of this new 
command strategy. 
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This work is divided into seven parts. In section 1, the introduction is presented. In section 2, the 
model of the DFIG is described. The DRWT has been discussed in section 3. In section 4, the description of 
the classical DPC technique is presented. Section 5 deals with the description of the DPC technique with the 
application of a genetic algorithm. Simulation studies are presented and discussed in section 6. We finished 
the paper with the conclusion. 


2. DFIG MODEL 
The mathematical model of the DFIG has been presented in detail in [31]-[33] by using the Park 
model. The equation of flux and voltages of rotor and stator DFIG are given as (1). 
d 
Vas = Rslas + a vas = WsWqs 
d 
Vs = Rslgs + a vas + WsWas 
d 
Var ma Rlar F a ar = Wr Wor 
d 
Vor = Rylgr + a var + @, War 


(1) 


The flux can be expressed as (2). 


Was = Mlar + Lslas 
Was = Mgr + Lslqs 
War = Lylar + Mlas 
Wor = Lylgr + M1q5 


(2) 


The reactive and active powers can be expressed as (3). 


3 
R= z Vaslas + Voslqs) 


3 (3) 
Q; = 5 Vaslas x Vaslqs) 


The torque is expressed as (4). 


dn 
Te=T, +] tf a (4) 


3. DUAL-ROTOR WIND TURBINE 

Traditionally, SRWT is a classical type of wind turbine used to this day in the production of 
electrical energy. The ideal maximum power coefficient of this type is 59%. So, the SRWT gives us a 
somewhart average coefficient. There is another type of wind turbine that gives a larger coefficient called 
DRWT. The latter has a coefficient estimated at 64%. Therefore, the DRWT improves the maximum power 
coefficient of 5% compared to SRWT [34]. It can be said that DRWT gives us more torque and mechanical 
power than SRWT. 

In the DRWT type, there are two turbines (Auxiliary turbine and main turbine). The block diagram 
of the DRWT show in Figure | [35]. The total aerodynamic torque of DRWT is the Auxiliary turbine plus the 
main turbine torque as shown by (5). 


Tprwr = Tr = Ty + Ta (5) 


Where: Tm: Main turbine torque. 
Ta: Auxiliary turbine torque. 
Tr: Total torque or DRWT torque. 
The aerodynamic torque of the auxiliary turbine are given by (6) [36]. 


ee! 
ad 3 
2143 


Ta .A. p. 1. R}. Cp.w2 (6) 


The aerodynamic torque of the main turbine are given by (7). 
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Ty = —.A. p.n. RS. Cp. we, (7) 


= 553 
243, 


With Ra, RM: blade radius of the main and auxiliary turbines, àa, Am: the tip speed ration of the main and 
auxiliary turbines, p: the air density and wm, wa the mechanical speed of the main and auxiliary rotors. 
Cp can be calculated as (8). 


1 0.035 
A+0.08B £341 








Cp (A,B) = (8) 
With ß is pitch angle 

The tip speed ratios for the main and auxiliary turbines are calculated through (9) and (10), 
respectively. 


_ Wa-Ra 
Ay = “a (9) 
pees (10) 
VM 


Where V; is the wind speed on an AWT and Vy is the speed of the unified wind on the main turbine. On the 
other hand, the essential element for calculating the tip speed ratio is wind speed on the main and auxiliary 
turbines. Obtaining the wind speed on the auxiliary turbine is straightforward. However, the calculation of 
wind speed on the main turbine requires further investigation. Based on (11), it is possible to estimate the 
amount of the wind speed at any point between the auxiliary and main blades. 


_ _ 1-VG=¢r) 2.x 
kena Sa + 2) (11) 


With x: the non-dimensional distance from the auxiliary rotor disk, Vx the velocity of the disturbed wind 
between rotors at point x and Cr the trust coefficient, which is taken to be 0.9. So, with respect to x=15, the 
value of the Vx close to the main rotor is computable (rotors are located 15 meters apart from each other) 
[34]. 


Ps_ref - 






Os ref —>» 
RSC 


Main Rotor 


Figure 1. Block diagram of DRWT with a DFIG. 


Figures 2-7 show the obtained simulation results. According to these figures, we notice that the 
mechanical powers extracted are adapted to the variation of the wind speed and the total mechanical power 
(Pm = 0.52 MW) at the moment t= 4 s is equal to the summation of the two mechanical powers secondary 
rotor (Pm2= 0.02 MW) and main rotor (Pm1= 0.50 MW). 
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Figure 2. Wind speed V1. Figure 3. Wind speed Vm. 
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4. CLASSICAL DPC METHOD 

DPC or direct power command using classical lookup table is the most used command strategy for 
DFIG-based DRWP systems. In this strategy, two hysteresis comparators are used to controlling the reactive 
and active powers. However, this method is a simple algorithm and reduced power oscillation, torque 
oscillation and harmonic distortion (THD) of stator current compared to field-oriented control (FOC). On the 
other hand, this method gives a fast response dynamic compared to the FOC strategy. Figure 8 shows the 
classical DPC method of DFIG driven by DRWT. 








Active, reactive power and sector estimation 


Figure 8. Classical DPC control. 
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The magnitude of rotor flux, which can be estimated by (12). 


2 -Ri Ë 
re Í (Vra r ia) (12) 


t š d 
9p = Í (V,g — R,- ip) £ 


The rotor flux amplitude is given by (13) and (14). 


$r = | Pra + prg (13) 


And 
p=% (14) 


The rotor flux angle is calculated by (15). 


6, = arctg (C2) (15) 


Reactive/active power is estimated using (16) and (17). 











P= Fa (bre) (16) 
Qs = — FE Org — bra) (17) 
Where: 

c=1- (18) 


In the classical DPC method, the level of reactive power hysteresis controller is two and the level of 
active power hysteresis controller is three. Figure 9 represents the reactive hysteresis controller, where this 
exists as the last output gives the last two values 0 and 1. On the other hand, the output of the active power 
hysteresis controller is 0, 1 and -1 as shown in Figure 10. The output of active and reactive power hysteresis 
comparators is the input of the switching table. The zone of the rotor flux is six-zone and the lookup table is 
shown in Table 1. 


Table 1. Switching table of DPC method 
N 1 2 3 4 5 6 

















Hq Hp 

1 6 1 2 3 4 5 

0 0 0 7 0 7 0 7 

-1 2 3 4 5 6 1 

1 5 6 1 2 3 4 

1 0 7 0 7 0 7 +0 

-l 3 4 5 6 1 2 

+APs/2 
-1 
Figure 9. Reactive power hysteresis comparator. Figure 10. Active power hysteresis comparator. 
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5. DPC-GA METHOD 

In this section, a novel DPC control was designed to regulate torque, active and reactive powers of 
the DFIG-based DRWT. This designed DPC strategy is based on pulse width modulation (PWM) and genetic 
algorithm. This proposed DPC method is a simple algorithm, fast response dynamic, robust control compared 
to classical DPC strategy and vector control. On the other hand, this proposed method reduced the reactive 
and active power of DFIG-based DRWT. So, the principle of the proposed strategy is a modification of the 
classical DPC method, where the classical switching table have been replaced by a PWM technique and the 
two hysteresis comparators of reactive and active powers has been replaced by to genetic algorithm. The 
proposed DPC strategy, which is designed to regulate the active and reactive powers of the DFIG-based 
DRWT is shown in Figure 11. 





(Gage eee mail 
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wind 





Active and reactive power estimation 


Figure 11. Block diagram DPC-GA of DFIG-based DRWT. 


6. SIMULATION RESULT 

The simulation results of the DPC-GA strategy of a 1.5 MW DFIG are compared with the 
conventional DPC strategy. The commands system was tested under deferent operating conditions such as a 
sudden change of load active and reactive powers. The performance analysis is done with torque, THD value 
of stator current, reactive/active powers. DFIG used for the simulations has the following parameters: P,=1.5 
MW, V,=398 V, f=50 Hz, f=0.0024 Nm/s, L;=0.0137 H, J=1000 Kg.m?, R=0.012 Q, R,=0.021 Q, L,=0.0136 
H, M=0.0135 H [37], [38]. Table 2 shows parameters of algorithm genetic. 


6.1. Reference tracking test (RTT) 

For the traditional DPC and designed DPC strategy, the reactive power, and active power track well 
their reference values (Psref and Qsref) as Shown in Figures 12 and 13. The active and reactive powers are 
decoupled from each other in the DPC-GA with a rapid time response, without overshoot, and with a 
minimal static error compared to the conventional DPC technique. Figure 14 shows the torque of the DPC 
and DPC-GA strategies. Figure 15 shows the stator current of the designed DPC strategy and the classical 
DPC. We can see that torque and current of the DFIG are proportional to the variation of reactive/active 
power reference values. Active power response comparing curves are shown in Figure 16. See figure the 
active power oscillations are significantly minimized when the DPC-GA strategy is in use. Figure 17 shows 
the zoom in the reactive power responses of both the DPC technique. It is found that the DPC-GA strategy 
exhibits smooth response and lesser oscillations in reactive power as compared to the conventional DPC 
technique. 

The DPC-GA technique reduced the torque oscillations compared to the classical DPC technique as 
shown in Figure 18. From the simulation results presented in Figures 19 and 20, it is apparent that the THD 
value of stator current for the DPC-GA technique is considerably reduced relative to the classical DPC 
technique and other strategies as shown in Table 3. Where, DPC-T2FLC is the direct power control with type 
2 fuzzy logic controller. DPC-MRAC is the direct power control with model reference adaptive control. FOC 
is the field-oriented control. DPC-NFC is the direct power control with neuro-fuzzy controller. 
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Table 2. Parameters of algorithm genetic 
Current iteration 51 
Number of variables 3 
Bounds lower [-10 -100 0] [200 2000 0] 
Population type Double vector 
Population Creation function Use constraint dependant defaut 
Population size Use default: 20 
Fitness scaling: Scaling function Rank 
Selection: Selection function Stochastic uniform 
Reproduction Elite count Use default: 2 
Cross over fraction Use default: 0.8 
Mutation: Mutation function Use constraint dependent defaut 
Crossover: crossover function Scattered 
Direction Forward 
Migration Fraction Use default: 0.2 
Interval Use default: 20 
Algorithm settings Intial penalty Use default: 10 
Penalty factor Use default: 100 
Hybrid function: hybrid function None 
Generations Use default: 100 
Time limit Use default: inf 
Stopping criteria Fitness limit Use default: -inf 
Stall generations Use default: 50 
Stall time limit Use default: inf 
Function tolerance Use default: 1e-6 
Nonlinear constraint tolerance Use default: 1e-6 
Output function: History to new window Interval: 1 
Display to command window Level of display: off 
User function evaluation Evaluate fitness and constraint functions: in serial 
Table 3. THD value (RTT) 
THD (%) 
DPC-PI 2.59 
Ref. [39] DPC-ANN 1.09 
Ref. [40] DPC-T2FLC 1.14 
DPC-NFC 0.78 
Ref. [41] FOC 3.7 
Ref. [42] SMC 3.05 
FSMC 2.85 
Ref. [43] DPC-MRAC 1.01 
Proposed strategy DPC 0.54 
DPC-GA 0.11 
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z zme | 
A Ps (DPC) 3 Qs (DPC) 
g 5 Ps (DPC-GA) a. Qs (DPC-GA) 
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Time (s) Time (s) 
Figure 12. Active power. Figure 13. Reactive power. 
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Figure 14. Torque. 
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Figure 17. Zoom in the reactive power 
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Figure 18. Zoom in the Torque 


6.2. Robustness test (RT) 
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Figure 19. THD value of current (DPC). 


Table 4. THD value (RT) 








THD (%) 
DPC DPC-GA 
Stator current 0.91 0.29 





In this part, the DFIG parameters have been intentionally changed such as the values of the 
resistances R, and R, are multiplied by 2 and the values of the inductances L, and L, are multiplied by 0.5. 
Simulation results are presented in Figures 21-26. As it’s shown by these figures, these variations present a 
clear effect on reactive power, torque, active power, and current curves and that the effect appears more 
important for the classical DPC technique than that with designed technique as shown in Figures 27-29. On 
the other hand, these results show that the THD value of current in the proposed strategy has been reduced 
significantly as shown in Figures 21-22. Table 4 shows the THD values of both techniques. Thus, it can be 
concluded that the designed strategy is more robust than the classical DPC technique. 
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Figure 20. THD value of current (DPC-GA). 


Figure 21. THD value of current (DPC) 
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In order to regulate the active and reactive powers of the DFIG-based DRWP system, a GA method 


is proposed. The results of simulations have shown that the application of a GA technique gives a good 
response of reactive and active powers. Power level oscillations are lower relative to classical DPC strategy, 
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which is reflected in the quality of the stator currents generated by the DFIG-DRWP. From the comparative 
study between the designed technique and the classical DPC strategy, it has been shown that the designed GA 
method is very effective in the stabilization of the system. Therefore, the designed strategy can contribute to 
expanding wind power utilization. 
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